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Dicke model predicts a quantum phase transition from normal to superradiant phases for a two- 
level atomic ensemble coupled with an optical cavity at zero temperature. In a recent pioneer 
experiment [Nature 464, 1301 (2010)], such a phase transition has been observed using a Bose- 
Einstein condensate (BEC) in an optical cavity. Compared with the original Dicke model, the 
experimental system features finite temperature and strong atom-photon nonlinear interaction in 
BEC. In this Letter, we develop a finite temperature theory for the Dicke phase transition of a 
BEC in an optical cavity, taking into account the atom-photon nonlinear interaction. In addition 
to explaining the experimentally observed transition from normal to superradiant phases at finite- 
temperature, we point it out that a new phase, the coexistence of normal and superradient phases, 
was also observed in the experiment. We show rich finite temperature phase diagrams existing in 
the experimental system by tuning various experimental parameters. We find that the specific heat 
of the BEC can serve as a powerful tool for probing various phases. 

PACS numbers: 37.30.+i, 42.50.Pq, 03.75.Hh 



The Dicke model, a textbook paradigm in quantum op- 
tics, describes a two- level atomic ensemble coupled with 
an optical cavity. It was first introduced to illustrate 
the importance of collective and coherent excitations for 
atoms induced by a single photon mode With in- 
creasing atom-cavity coupling, the Dicke model predicts 
a quantum phase transition from a normal phase (NP) 
to a superradiant phase (SP), where both the atomic en- 
semble and photon acquire macroscopically collective ex- 
citations [2|-|4[ . However, due to the 'no-go theorem' orig- 
inating from the Thomas-Reiche-Kuhn sum rule for the 
oscillator strength, this phase transition cannot be re- 
alized in typical cavity quantum electrodynamics [H, [6[ . 
The experimental breakthrough for observing the quan- 
tum phase transition only occurs recently using the mo- 
mentum eigenstates of a Bose-Einstein condensate (BEC) 
coupled with an optical cavity @, Q • I n this pioneer ex- 
periment, there exists an nonlinear atom-photon interac- 
tion which, induced by the optical lattice potential 
and greatly enhanced by the large atom number TV, can 
reach the same order of the effective cavity frequency 
and even go beyond. In this strong interaction regime, 
rich dynamical properties Io|, Ell as well as new quan- 
tum phase transitions [12| have been predicted at zero 
temperature. 

Although such ideal quantum phase transitions should 
occur at absolute zero temperature [l3[ in principle, re- 
alistic experiments used to observe the quantum phase 
transition must be performed at finite temperature. For 
instance, the typical temperature of the BEC in 0,0 for 
observing the Dicke phase transition is ~ 50 nK. At finite 
temperature, thermal fluctuations may induce new exotic 
phenomena beyond the prediction of the zero tempera- 



ture theory [14[ . A well-known example is the lost of long 
range superfluid order in low dimensions (two or one) 
at any temperature [15], where the superfluid physics 
is characterized by the Berezinskii-Kosterlitz-Thouless 
transition [l6|, [ItJ- Therefore it is crucially important 
to investigate the Dicke phase transition at finite tem- 
perature to fully understand the realistic experiment. 

In this Letter, we present a finite-temperature field the- 
ory for the Dicke phase transition in a BEC coupled with 
an optical cavity to understand the recent breakthrough 
experiment in this system. Our main findings are the 
following: 

(I) The experimentally observed phase diagram for the 
normal-superradiant phase transition, (i.e., Fig. (5) in 
Ref. l 7]) can be well understood in our finite-temperature 
theory. More interestingly, we show that a new phase, 
called CEsp, can be identified in the experimental phase 
diagram in [7[. In the CEsp phase, NP and SP coexist, 
but the SP is stable and the NP is met ast able. 

(II) By varying various physical parameters (temper- 
ature, the coupling strength, the nonlinear interaction, 
etc.), we show there exist rich finite-temperature phases 
in the experimental system, including NP, SP, CEsp, 
CEnp (coexistence of stable NP and metastable SP), as 
well as dynamical unstable phases (US). In certain pa- 
rameter region, there is a four-phase coexistence point. 

(III) We show, both analytically and numerically, that 
the specific heat of the BEC increases rapidly in the 
SP / CEsp (exponential increase at low temperature) with 
increasing temperature, and has a large jump at the criti- 
cal transition temperature where the system becomes the 
NP/CEnp- Therefore the specific heat of the BEC may 
serve as a powerful tool for detecting different phases in 
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FIG. 1: (Color online) (a) A BEC with the 87 Rb atoms in- 
teracts strongly with a high-finesse optical cavity, (b) A 
four-level model is considered by introducing the zero mo- 
mentum state |0, 0) (= \p x ,Pz) with p x and p z being the mo- 
menta in the x and z directions, hereafter), the excited states 
|d=fc, 0) and |0, =b/c) , and the symmetric superposition of states 
|=b/c, dz/c). In the dispersive limit, the excited-state levels can 
be eliminated adiabatically, and an effective two-level system 
with the zero momentum state and the symmetric superposi- 
tion of states is thus formed. 

the temperature-driven Dicke phase transition. 

Figure Q] shows the experimental setup in which all 
atoms in a 87 Rb BEC interact identically with a single- 
mode photon induced by a high- finesse optical cavity [7j . 
In the experimental scheme, the atom momenta are used 
to define the spin states \\) = |±fc, ±fc), \\) = |0,0) with 
the corresponding SU(2) collective spin operators S+ = 

Sl = It)« Ul and S z = Eidt)^ (t| - \i) u (II), as 
shown in Fig. [Tfb). As a result, not only the 'no-go 
theorem' is overcome, but also the superradiant-normal 
phase transition condition can be satisfied (191421). since 
the smaller energy scale of the effective two levels can be 
achieved. Under the new spin basis, the dynamics of the 
atom-cavity system are governed by the Hamiltonian [?} 

h = (co + ^smH + Y Sz + 7W {lp + ^ ){s+ + s - ] - 

(i) 

Henceforth we set ft = 1. ^ denotes the photon cre- 
ation operator. The effective cavity frequency uo = 
— A c + NUo(l + S)/2 with A c = uj p — uo Cl where uo c is 
the cavity frequency, uo p is the pump laser frequency, 
H = 3/4, and Uq = g\j{uj v — uj a ) with g being the 
coupling strength between a single atom and the pho- 
ton and uj a being the atomic transition frequency. The 
nonlinear atom-photon interaction U = NEUo, which is 
greatly enhanced by the large atom number TV, and can 
reach the same order as the effective cavity frequency uo or 
even beyond. Therefore a strong nonlinear atom-photon 
interaction regime is accessible in experiments. The ef- 
fective atomic frequency ujq = 2uo r with the atomic recoil 
energy uo r = k 2 /2m. The collective coupling strength 
g = g Qy/N /2(uj p — w a ), and Q is the maximum pump 
Rabi frequency. 



The finite-temperature properties of the BEC-cavity 
system can be obtained by calculating the partition func- 
tion of the system through the imaginary-time functional 
path-integral approach. In this procedure, we rewrite 
the collective spin operators in the Hamiltonian (pQ) as 

S z = T,iLi( a i a i - llli) and S + = J2^=i a hi usin & 
Fermi operators d\(cxi) and which obey the anti- 

commutation relation {a\,aj} = {jj^Jj} = After a 
straightforward calculation, the effective partition func- 
tion can be written as [l8[ 

Z_ = J[d V ]exp(-S) 

Z J[drj}eM-SfY [ ' 

where [drj] is the functional measure. In 
Eq. (J2j), the free action of the bosonic field 
Sf = Jo dr Eti K* (r) d T a, (r) + 7* (r) d rll (r)] + 
Jq dri/j* (r) (d T + uj) ip (r), where r = it, d T = d/dr, 
(3 = 1/ {ksT) with ks being the Boltzmann constant 
and T being the system's temperature. The total action 
S is expressed, in the basis of <&i (r) = [7^ (r) , (r)] T , 
as 

S = So (V, r) + f dr*t (r) G ^*) % (r) , (3) 
Jo 



where So (-0, 0*) = J drip* (r) <9 r (r) and 
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FIG. 2: (Color online) The scaled mean-photon number as the 
functions of the detuning A c and the square of the collective 
coupling strength g for the temperature T = 0.00 nk (a), 
T = 0.13 nk (b) T = 50.00 nk (c), and T = 150.00 nk (d). 
The nonlinear interaction U — 1.70 x 10 3 cj r and the recoil 
frequency uj r = 2tt x 3.71 kHz. 1 nk +± 130.00 Hz. 



3 



005 0.1 2 4 6 8 




c 


) 2 


4 


C 


) 5 


10 


1 ' 




















0.8 








US 


SP 


- 0.6 

S 

H 0.4 


NP 


CE SP 








0.2 

rv 


(a) 






NP 

(b) 


CE sp 



25 



30 

U/co 



3525 



30 

U/co 



35 



FIG. 3: (Color online) The scaled mean-photon number as the 
functions of the temperature T and the collective coupling 
strength g for U = 10.00 cj r , A c = -20.00o; r (a) and U = 
31.10 cj r , A c = 20.00^ r with uj r = 2tt x 3.71 kHz. 



FIG. 4: (Color online) The scaled mean-photon number as the 
functions of the temperature T and the nonlinear interaction 
U for g = 3.00 cj r , A c = 20.00cj r (a) and g = 6.00 cj r , A c = 
20.00cj r with uj r = 2tt x 3.71 kHz. 



We first integrate out the Fermi field &i (r) (thus the 
atom spin degrees of freedom) in Eq. (j2j), and then use 
the standard stationary phase approximation [22], i.e., 
SSfy*,^) /6ip(r) = and SS (ip*, ip) /Si/j* (r) = 0, to 
obtain the finite temperature equilibrium phase diagram 
in the large atom number, which is governed by the fol- 
lowing equation 

ujifj = F(^*,^)tanh[/3G(^*,^)]^, (5) 

and its complex conjugate. Here F(tp*,il;) = 
[(f + U + 4 5 2 ] /[2 V / (^f + ^^) 2 +4^^] 

and G(V*,V) = ^(f + ^V'*V') 2 +4f^V/2. It is 
straightforward to find that Eq. (j5j) has a trivial solution 
tjj = ^* = and the non-trivial solution ip (r) = 
^* (r) = ±^0? satisfying a nonlinear equation 

(^|p7W = tanh (^)' (6) 



where ( = y (w + Ui/>$) + 16# 2 ^o, and = 
{ipty) /N is the scaled mean photon number. If ip = 
ijj* =0, the system is in the NP with no collective 
excitation. However, the existence of non-trivial solu- 
tions shows that the macroscopic collective excitation 
for both atoms and photon occurs and thus the sys- 
tem locates at the SP. Finally, the stable conditions 
5 2 S{^\^)/5^ 2 > and S 2 S (^*, ip) /S(^) 2 > de- 
termine the real solution of Eq. (j5j) and thus the phases 
of the BEC-cavity system at finite temperature. 

When U = 0, Eq. ® reduces to 

uy/ul + Wg 2 ^ 2 /4g 2 = tanh (/3 y/u* + 1 6g 2 ^ 2 /^) , 
the result for the original Dicke model. In this case, only 
SP and NP exist, even at finite temperature [l8[. When 
the temperature T — ^ 0, tanh (/3C/4) = 1, and Eq. © 
becomes 2cj£/[(co>o + U^q) U + 8g 2 ] = 1, from which rich 



zero-temperature phase diagrams have been obtained 
With increasing temperature T, the function 
tanh(/?£/4) decreases. As a consequence, Eq. (|6]) can 
possess some new physical solutions of t/Sqj leading to 
new phases at finite temperature. 

We solve the nonlinear equation (j5j), together 
with stable conditions 5 2 S (tjj* , ip) / Stjj 2 > and 
5 2 S (^*, iji) /^(^*) 2 > 0, for various different physical pa- 
rameters and calculate the scaled mean photon number 
(ipty) /N = t/Sq. In Fig. 2, we plot the scaled mean- 
photon number i/jq as the functions of the detuning A c 
and the square of the collective coupling strength g for 
different temperatures (a) T = 0.00 nK, (b) T = 0.13 nK, 
(c) T = 50.00 nK, and (d) T = 150.00 nK. To compare 
with the experiment, we also choose the other parameters 
to be the same as those in Ref. 0, that is, the nonlinear 
interaction U = 1.70 x 10 3 cj r , and the recoil frequency 
uj r = 2tt x 3.71 kHz. When T = 0, due to the existence 
of the strong nonlinear interaction, the CEnp phase is 
found, apart from the NP and SP, as shown in Fig. 2 
(a). In the CEnp phase, ~ 0, and almost no pho- 
ton signature is generated. However, at the experimen- 
tal operating temperature T = 50.00 nk, numerical solu- 
tion of Eq. (|5]) shows that the finite-temperature critical 
point that separates the NP and SP agrees with that in 
the experiment. More interestingly, the CEnp phase be- 
comes the CEsp phase, in which a weak photon signature 
(compared with the strong superradiant regime) emerges. 
We note that such weak photon signature for the CEsp 
phase has already be observed in the experiment in the 
same parameter region, although it has not been explic- 
itly pointed out. The transition temperature from the 
CEnp phase at zero temperature to the CEsp at 50 nk 
is T p = 0.13 nk, as shown in Fig. 2(b). With increas- 
ing temperature T (=150.00 nk), the SP region becomes 
narrower due to the suppression of the superradiance 
by thermal fluctuations, as expected. However, a new 
phase diagram including the NP, SP, CEsp, CEnp and 
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FIG. 5: (Color online) The specific heat C v as a function of the 
temperature T for different nonlinear interaction U, detuning 
A c and collective coupling strength g. Insert: comparsion be- 
tween the numerical simulation (NS) and the analytical result 
(AR) in Eq. J7j) for the specific heat C v . 



US (5 2 S(^\^) /Sip 2 < and S 2 S (^*, ip) /S(ip*) 2 < 0) 
is observed, as shown in Fig. 2(d). 

In a realistic experiment, the experimental parameters 
such as the detuning A c , the collective coupling strength 
g, the nonlinear interaction U as well as the temperature 
T can be controlled independently. For example, A c and 
g can be tuned by adjusting the frequency and power of 
the pump laser, respectively. The nonlinear interaction 
can be controlled by tuning both the frequency of the 
pump laser and the atom number. The temperature can 
be manipulated by controlling the evaporative cooling 
process for the BEC [23[. In Fig. 3, we plot the scaled 
mean-photon number t/Sq as a function of the temperature 
T and the collective coupling strength g. While in Fig. 
4, we plot the scaled mean-photon number t/Sq as a func- 
tion of the temperature T and the nonlinear interaction 
U. These figures show that rich finite- temperature phase 
diagrams exist in different parameter regions, induced 
by the competition among these parameters. Interest- 
ingly, we find a four-phase coexistence point in certain 
parameter region, as shown in Fig. Hfb). It should be 
emphasized again that these exotic phases arise from the 
interesting nonlinear interaction U and the finite tem- 
perature. If U = 0, only the SP and NP can be found 

m 

These rich finite temperature phases and the phase 
transitions can be detected in experiments by measur- 
ing not only the mean photon number, but also some 
the thermodynamic quantities, such as the specific heat 
per atom Cy = N j^ bT i -§m (In Z) in the BEC. In the 
NP, Cy can be obtained analytically, yielding Cy P = 
ujI [l - tanh 2 (/3cj /4)] /(Sk B T 2 ) [2], which reaches the 
maximum at certain temperature (the red dashed-dotted 
line in Fig. 5). However, the explicit expression for Cy 
cannot be obtained in the SP. Nevertheless, in the region 
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where Co is the value of ( = y (ojo + Utp 2 ,) + 16 g 2 ^ 
at T = 0. This expression shows that the specific heat 
Cy P in the SP increases rapidly (exponential increase at 
low temperature), which agrees well with the numerical 
results (see the insert of Fig. 5). At the critical temper- 
ature T c , at which the system enters the NP, the specific 
heat Cy has a large jump. This step behavior is quite dif- 
ferent from the behavior of the scaled mean photon num- 
ber ?/5q, which varies smoothly when crossing the critical 
point. Therefore the temperature-driven phase transition 
from SP to NP can be detected using the specific heat 
Cy. In addition, the behavior of specific heat Cy in the 
CE S p/CE NP is similar to that in the SP/NP. Moreover, 
for the transition from the CEsp to CE^p phases, such 
a large jump of Cy still exists, as shown by the green 
dotted line of Fig. 5. 

In summary, we develop a finite temperature theory for 
the Dicke phase transition of a BEC in an optical cavity. 
Our theory not only explains well the recent experimen- 
tal observation of the phase transition from the SP to 
NP, but predicts rich new phases that either have been 
observed (but not pointed out) or may be observable by 
tuning experimental parameters. Our study is crucially 
important for understanding the Dicke phase transition 
in a realistic experimental system and may have signifi- 
cant application in quantum computation, quantum op- 
tics, etc. 
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